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Abstract 

The Coulomb pseudo-potential /U* is extracted by fitting the numerically calculated transition 
temperature Tc of the Eliashberg-Nambu equation which is extended to incorporate the narrow- 
band effects, that is, the vertex correction and the frequency dependence of the screened Coulomb 
interaction. It is shown that even for narrow-band superconductors, where the fermi energy ei;' is 
comparable with the phonon frequency ujph, the Coulomb pseudo-potential is a pertinent parameter, 
and is still given by /i* = + nln{eF/ujph)], provided tOph is appropriately scaled. 
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It is generally agreed that the superconducting properties of the fulleride superconductors 
can be understood in terms of phonon-mediated s-wave pairing. The relevant physical pa- 
rameters are: the logarithmically averaged phonon frequency uoin ~ 0.1 eV, the Fermi energy 

~ 0.2 — 0.3 eV, the dimensionless electron-phonon coupling constant A ~ 0.7—0.8, and the 
Coulomb repulsion /i = Vciyj = 0)Np ^ 0.3 — 0.4, where Vc{uj) and Np are, respectively, the 
screened Coulomb repulsion and electronic density of states at the Fermi level 0] . There- 
fore, in fuUerides we have ep ~ (^in unlike the conventional metals where ep ^ uin- In this 
case, the Coulomb pseudo-potential /i*, which is given by fi* = + fi\n{ep/ujph)], is not 
much reduced from ft as ep becomes comparable with the average phonon frequency Uph (see 
below for further discussion on Uph and Win). The observations that the value of \n{ep/u!ph) 
is large, typically 5 — 10 in conventional wide-band superconductors, and that the frequency 
dependence of the screened Coulomb interaction is negligible in the phonon frequency region 
are crucial in introducing the concept of /i* 0. Both of these observations do not hold in 
the fullerenes because ~ uin, and, consequently, it is not clear whether the concept of 
Coulomb pseudo-potential is still valid in this case. In the present Letter, therefore, we study 
the Eliashberg equation extended to narrow-band superconductors to check if the Coulomb 
pseudo-potential remains a pertinent parameter as in the wide-band superconductors, that 
is, if Tc is determined solely by uin, A and /i*, where /x* is independent of A and contains all 
the effects of the Coulomb repulsion. We first obtain the superconducting transition tem- 
perature Tc by numerically solving the Eliashberg-Nambu equation. Then, the calculated Tc 
values were fitted with a modified Tc form which is a function of A and /i, in order to extract 
fi*. We find that even for narrow-band superconductors the Coulomb pseudo-potential is a 
physically relevant concept, and is still given by /i* = /i/[l -|- filn{ep/u!ph)], provided that 
Ljph is appropriately scaled. 

For wide-band superconductors, the frequency dependence of the screened Coulomb in- 
teraction is not important. The reason is that Vc{uj) does not change much on the phonon 
frequency scale of Uph when ep ^ Uph, since the frequency scale over which Vc{uj) varies is ei? 
and the pairing kernel becomes negligible beyond the region of width Uph around the Fermi 
surface. This does not hold for narrow-band superconductors, and it may be important 
to take uj dependence of Vc{uj) into account. We do this by modeling the Coulomb repul- 
sion in terms of the onsite Hubbard repulsion and calculating its screening self-consistently. 
Another complication for narrow-band superconductors is that the Migdal theorem does 
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not hold and the vertex correction may become important. The vertex correction is in- 
cluded in the present work following the Nambu method [^. These extensions, that is, 
(a) frequency dependence of the screened Coulomb interaction and (b) vertex correction, 
are included in the Eliashberg equation. Gunnarsson and his coworkers pointed out that 
the procedures of screening and renormalizing away the high energy states beyond a cutoff 
frequency do not commute in narrow-band superconductors and there arises a subtle 
problem in the concept of Coulomb pseudo-potential. This problem does not occur in the 
present approach. Including the full frequency dependence of the self-consistently screened 
Vc{uj) in the Eliashberg equation unnecessitates the introduction of the cutoff frequency 
in manipulating Vc{uj) in the Eliashberg Tc equation 0, which becomes ill-defined in the 
narrow-band superconductors. 

It is straightforward to derive the Eliashberg-Nambu (EN) equation for finite bandwidth 

superconductors in the Matsubara frequency 0: 
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l^H I '^ ={X-{vi(p^ + V2W^) + \+V:iW^], (1) 

P m JWl + 61 



where pn = 27iT{n + 1/2) is the Matsubara frequency, 9n = ta.n^^ epj {^Zn\Jpn + ^n) 
and T is the temperature. Wn = PnZn and (pn = ^nZn, where A„ = A{ipn) and Zn = 
Z{ipn) are, respectively, the superconducting pairing and the renormalization functions when 
analytically continued into the real frequency {ipn uj + i6) , and 

A+ = \ph{n - m) - Xch{n - m) + \sp{n - m), 

A_ = \ph{n -m) - Xch{n - m) - A^p(n - m). (2) 

Here, \h{n — iti) is the pairing kernel due to the electron-phonon interaction at the Mat- 
subara frequency {pn — Pm) = 2ti {n — m)T given by 

where Q!^F(f2) is the phonon spectral function. Xch and Xgp are, respectively, the interactions 
in the charge and spin channels due to the Hubbard repulsion U, and are determined self- 
consistently as 

Xch{n-m) = UNF{l/2-{Xn + Xs) + iXn + XsfH'^ + l/iXn + Xs)]}, 

X,p{n - m) = UNf{1/2 + {xn - Xs) + {Xn - XsfH^ - l/{Xn - X.)]}, (4) 



where Xnin — m) and XsiP' — f^) are the dimensionless susceptibihties from, respectively, the 
normal and pairing processes given by [|^ 

Xn{k) = — -r}^ 



The finite bandwidth is explicitly incorporated through the factor 6, which is equal to 7r/2 
in the conventional infinite bandwidth case. Finally, the vertex corrections are included by 
extending the Nambu method M. The vertex functions are given by 
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where ujpi is the plasma frequency which is taken to be 1 eV. Without the vertex correction, 
which corresponds to f i = 1 and f2 = "^3 = 0, Eq. (|l|) is of the same form as in the theory 
used to study spin fiuctuation effects on superconductivity p. Self-consistent solution of 
Eq. together with Eqs. (H)-® yields Z{ip„) and A(zp„) in the Matsubara frequency. 
Details of the formulation will be reported elsewhere 0. 

It is simple to calculate from the EN equation of Eq. ([l|). It is identified as the highest 
temperature at which the pairing function A is non-vanishing. The input parameters for 
solving the coupled EN equations are: a^F(fi) (which determines uin and A), U, and ep- 
Note that the long range part of the Coulomb interaction Vq is not included in the present 
study. This may be justified for the case where the average of Vg is small compared with 
the onsite Hubbard repulsion U: (yq)avel^ ~ Anne'^/kpU ^ 1. Alternatively, U may be 
considered as a partially screened value due to the metallic screening of the long range part 
of the Coulomb interaction. The uoin is given by 

'2 /-"^ „ 1 
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where A = NpV, /i = UNp, and ^l/ = 2 /o°° dna^F{9)/n. We take a^F{9) as a sum of three 
truncated Lorentzians appropriate for fuUerene superconductors 0. This choice of a^Fip) 
gives the logarithmically averaged phonon frequency ujin = 0.094 eV. is varied between 



and 0.4 and A between 0.4 and 1.2 to calculate Tc, from which to extract an extended Tc 
formula and the Coulomb pseudo-potential fi*. 

The EN equation was solved by iterations, for various values of A and /i to find T^. Self- 
consistency is achieved with a few tens of iterations except when T is close to T^. These 
results are used to extract a modified formula following the same procedure as McMillan 
and Allen (MA) fTg. The MA equation is given by 



ir = — exp 
1.2 ^ 



1.04(1 + A) 



(8) 



A-/i*(l + 0.62A)_ 

which was extracted by fitting to numerical calculations of the Eliashberg equation. The 
McMillan-like procedure for extracting formula proceeds as follows JTOf: Firstly, for the 
case of zero Coulomb potential, /i = /i* = 0, we estimate the A dependence of Tc by taking 



In 



b c 
"+A + P' 



(9) 



where a, b, and c are the constants to be determined by fitting to the exact numerical re- 
sults obtained by solving EN equation of Eq. (|l]). The MA case corresponds to a, 6 7^ 
and c = 0. Non-zero c introduces higher power of 1/A in the formula due to the narrow 
bandwidth effects considered here, that is, the vertex correction and frequency dependence 
of the effective Coulomb interaction. McMillan used an approximate solution of the Eliash- 
berg equation to guide the choice of the functional form of Eq. (|^). From a solution of a 
version of a square-well potential model, the form of In{ujin/Tc) = a + b/X was anticipated 
||lO| . The simplest correction coming from beyond the simple model will be of the form 1/A^ 
as incorporated in Eq. (^. For the parameter values appropriate for the fuUeride supercon- 
ductors, we obtain a = 1.13, b = 0.765, and c = 0.033. This should be compared with the 
MA case, a = 1.22, b = 1.04, and c = 0. The obtained values for a, b, and c depend weakly 
on the particular choice of ep/^in and a^F. 

Secondly, to include the effects of non-zero Coulomb interaction on (/i 7^ 0), we take 
the form 

aA^ + b\ + c 
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A2 - ij,**{ij,,X)' 

as with MA. Here, a, b, and c have already been determined from the first fitting process, 
and jj,** depends on and A. For the MA case, n** = /i*(l + 0.62A), where 



H* = fj./[l + fj.ln{eF/ujph)] 



(11) 
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FIG. 1: Determining S = u)ph/u)in by requiring be independent of ^, when the bandwidth 

B is 0.5 eV. A is varied from 0.4 to 1.2, and, for each A, four /i values, 0.1, 0.2, 0.3, and 0.4, were 
taken to calculate fi** /fi*. (a) When S = 1. It is clear that for each A, the four calculated points 
corresponding to four different fi values are scattered in the plot, (b) When S = 2, the scattered 
4 points in the figure (a) scale to collapse on to a single point for each A. 



as given by Morel and Anderson [11|. We put 
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where /(A) is a function of A to be determined by fitting, and fi* is assumed, in the present 
case as well, to take the same Coulomb pseudo-potential form of Eq. ([Tl|). Then, /(A) can 
be written as 
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FIG. 2: S = ojph/uji n as a function of the bandwidth B. S was determined by requiring ji** / fi* be 
independent of ^ for a given bandwidth. The + and x are, respectively, calculated values for S 
and epjujph- The solid and dashed lines are a guide to eyes. 

We then write 

CUph = SuJin, (14) 

where 5* is to be determined by requiring that ii** / ii* = /(A) be independent of fi. We 
calculate, for a fixed A, n** / n* for several values of /i. fi**/^*, in general, will depend on /i 
and be scattered in the plot as shown in Fig. 1(a) for S = 1. For an appropriate value of 
S, however, the data will scale to collapse on to a single point for each A. This is found to 
occur at S* = 2 when ep = 0.25 eV as can be seen in Fig. 1(b). Fitting the curve of Fig. 1(b) 
yields /(A) = = O.eiSA^ + 0.244A + 0.244. 

This scaling behavior gives reliability to our procedures. Consequently, Tc can be written 

as 

aX^ + bX + c 



MA case corresponds to c = c' = 0. The concept of Coulomb pseudo-potential is found to 
be valid and important even for a narrow-band superconductor of ep/ojph ~ 1 where retar- 
dation effects are much weaker compared with the conventional wide-band superconductors. 
We should note here that although we determined 5* = Uph/ojin by requiring fi**/fi* be in- 
dependent of yU, Uph should not be considered an additional fitting parameter. S remains 
the same for a given bandwidth irrespective of A and /x, and gives the appropriate Coulomb 
pseudo-potential value fi* for given /x and bandwidth. 



We then vary the bandwidth B = 2eF to find how S varies accordingly. For each B, 
we repeated the fitting procedures described above to determine S = cuph/ujin which gives 
the best scahng behavior. as a function of B determined this way is shown in Fig. 2. S 
increases as the bandwidth is increased, and then saturates around the value of 4 when B 
is larger than 1 eV for uin = 0.094 eV. Note that in the Eliashberg theory, the Coulomb 
pseudo-potential /i* of Eq. ( ]TT| ) is often rescaled with a cutoff frequency c<Jc ~ 5 — 10 ud, 
where is the Debye frequency. Uc corresponds to Uph in the present work. The value 
of S" ~ 4 for i? > 1 is consistent with the uj^/ujo ~ 5 often adapted for conventional wide- 
band superconductors 0. Increase of 5 as -B is increased does not mean the less effective 
retardation, because the increase of S is slower than that ofep- This is shown with the crosses 
in the plot, ep/ujph increases as B is increased, and the retardation, which is determined by 
^p/^ph as given by Eq. (|ll|), is more effective for wide-band superconductors. 

To summarize, we investigated whether the Coulomb pseudo-potential /i* is still a perti- 
nent concept in the narrow-band superconductors. Because the narrow bandwidth implies a 
negligible retardation effects (/i* ~ yu) , it is not clear whether the /i* remains a relevant and 
useful parameter for the narrow-band superconductors like the fullerenes. To answer this 
question, we solved the extended Eliashberg equation which incorporates the finite band- 
width effects, namely, the frequency dependence of the effective interaction between electrons 
and the vertex corrections, and obtained the transition temperature as the dimensionless 
electron-phonon coupling constant A and the Coulomb repulsion ^ are varied. The obtained 
numerical results were used to find a Tc formula which is valid for narrow-band superconduc- 
tors and to extract the Coulomb pseudo-potential yU*. It was found that Tc is given by the 
form of Eq. ( pTSf ) where /i* contains all the effects of the Coulomb repulsion and is indepen- 
dent of A. Therefore, the concept of Coulomb pseudo-potential continues to be valid even for 
narrow-band superconductors. Another interesting and important outcome of the present 
study is that the Coulomb pseudo-potential is still given by /i* = + n\n.{eF/ojph)] as in 
the conventional wide-band superconductors, provided ojph = Suin is appropriately scaled. 
Scaling of Uph and /i* as Ci;' is varied is also computed. Because we included the vertex 
correction and frequency dependent Coulomb interaction in the present work, there is no a 
priori reason that the Coulomb pseudo-potential take the same form. We established that 
the Coulomb pseudo-potential continues to be a pertinent parameter and is still given by 
the Morel-Anderson form in the narrow-band superconductors. 
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